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Abstract 

In the hadroproduction of charm (or heavy flavours in general) in the con- 
text of string fragmentation, the pull of a beam remnant at the other end 
of a string may give a charm hadron more energy than the perturbatively 
produced one. The collapse of a low-mass string to a single hadron is the 
extreme case in this direction, and gives rise to asymmetries between lead- 
ing and non-leading charm hadrons. We study these phenomena within 
the Lund string fragmentation model and improve the modelling in part 
by a consideration of hadroproduction data. Applications include heavy 
quark production in any collision between hadron-like particles such as 7p 
at HERA and pp at HERA-B or the LHC. 
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1 Introduction 



Several experiments at fixed-target energies have observed asymmetries between leading 
and non-leading charmed hadrons in 7r~p collisions M. In H we study this effect in the 
context of the Lund string fragmentation model || and find that it is possible to improve 
the model to obtain a good description of the asymmetry data. In this model, collapses 
of light colour singlets mainly account for the asymmetry. 

We also study single charm spectra and particle/anti-particle correlations and find 
good agreement between the model and single charm spectra, but fail to describe the cor- 
relation data from E791 [|J. This is in contrast to NLO calculations ||, which adequately 
describe single charm spectra for non-leading particles and some correlations, but fail to 
describe the asymmetry and leading-particle single charm spectra. 

A related effect included in the string model is the non-perturbative 'beam drag'. In 
this model the produced heavy quarks are connected to the beam remnants by colour 
singlet strings. Through soft interactions with the beam remnant the charm quark can 
gain energy and momentum to produce a heavy hadron at a larger rapidity than that of 
the heavy quark. 



2 Model aspects 

The asymmetry between D + and D~ is described by the asymmetry parameter 

s o-(D-) - ct(D+) 

which is observed to be an increasing function of xp in 7r~p collisions. In fig. [l] experimental 
results are compared to a modified version of the string model as implemented in the 
Pythia 6.1 § event generator. In the Lund string fragmentation model the produced 
charm quarks are connected to the beam remnants by colour singlet strings as e.g. in 
fig. 0. The colour of the incoming anti-up quark is inherited by the outgoing anti-charm 
quark, which forms a colour singlet together with the pion beam remnant down quark. 
The distribution of c-d string masses near threshold is shown in fig. [| both using the 
default parameters and with some modifications. 

Depending on the mass of a string, different hadronization mechanisms are used: 

Large mass (M string > m c + m q + 1 GeV) Ordinary string fragmentation with a contin- 
uum of phase-space states. 

Low mass The colour-singlet system forms a cluster that undergoes two-body decay into 
two hadrons. 

Very low mass If the invariant mass of a cluster is so small that two-body decay is 
not possible, the cluster is forced to collapse into a single hadron with compatible 
flavour contents. This gives rise to an asymmetry favouring leading particles. 

The qualitative nature of the asymmetry is thus described by the string model. How- 
ever, the quantitative behaviour depends on the parameters of the model. The model 
aspects considered are described in more detail in ||, so here we just give a brief sum- 
mary: 
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Figure 1: Experimental results on asymmetries compared to the modified version of 
PYTHIA 6.1 described in this note. Data is from 7r~p fixed-target experiments, where 
x-p > is the pion fragmentation region. 




Figure 2: String configuration in a ir p collision with quark fusion charm production. 




Figure 3: Distribution of c-d colour singlet masses in a 7r~p collision using different 
parameters of the model. The dashed curve in the first figure is the distribution of 
invariant mass of the hadronic system produced from the colour singlet string or cluster. 
Below the two-particle threshold only the one-particle states D~ and D*~ are available. 
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beam remnant 

*-ud 

beam remnant 

Figure 4: The ud diquark of the proton remnant has, at least in the average, a larger 
energy /momentum than the lone u quark. Because the diquark is in a colour singlet with 
the c-quark, in average, we will have (y D ) > (y c ) and (y D ) > for y be amr e mnant > 0. 

• The masses used for the light and heavy quarks regulate the fraction of collapses. 
More collapses gives rise to a larger asymmetry. We use m c + m q = 1.5 GeV + 
0.33 GeV, i.e. constituent masses and the charm mass that gives the best fit to the 
total charm cross section. 

• Width of the primordial k± distribution. A larger intrinsic k± of the partons en- 
tering the hard interaction gives rise to an asymmetry also at larger p±_ which is in 
agreement with data. We use a Gaussian with (k±) = 1 GeV. 

• Beam remnant distribution functions. The energy-momentum of the beam remnant 
has to be split between the constituent quarks. This aspect is varied by using 
different beam remnant distribution functions parameterizing the fraction of energy- 
momentum taken by the constituents. We use an intermediate scenario. 

• Threshold behaviour between collapse and decay of a cluster. The transition be- 
tween cluster collapse and cluster decay as the cluster mass is increased can be slow 
or fast. We use an intermediate scenario. 

• Energy conservation in collapse. When a string with a very low mass collapses into 
a hadron, both energy and momentum cannot be conserved simultaneously. We 
consider different ways to shuffle energy between strings and partons in an event, 
but find that observable quantities are not sensitive to this aspect. 

3 Beam drag 

Beam remnant drag is a different, less obvious, source of asymmetry. In an independent 
fragmentation scheme the rapidity of a charm quark is preserved, on the average, during 
fragmentation. This is not necessarily so in string fragmentation, where the produced 
hadron inherits energy and momentum from both endpoints of the string. The resulting 
hadron will then have a rapidity that, on the average, is shifted in the direction of the 
other end of the string. An example with a proton beam remnant is shown in fig. f|. 

This is why the asymmetry for negative x-p in fig. [I] is negative, even though D~ is a 
leading particle also in this region. This behaviour is strongly dependent on the choice of 
beam remnant distribution functions, see 0. 
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Figure 5: B° (dashed) and B (full) distribution and the resulting asymmetry, A = 
(a(B°) - a(B°))/(a(B°) + tr(B°)) at HERA-B (top) and the LHC (bottom). 

4 Applications 

We have tuned the parameters of the model from a consideration of fixed target charm 
production data. Other applications include any collision between hadron-like particles 
e.g. 7p, pp or pp. This sections gives just a few examples. 

Asymmetries in the B°B° system at pp colliders must be considered as a background 
to asymmetries from CP- violation 0. Fig. |^ shows the bottom asymmetry as a function 
of rapidity for HERA-B and LHC energies (pp at 40 GeV and 14 TeV respectively). At 
HERA-B the effect is as large as in fixed-target charm experiments and a measurement 
of the effect would be feasible. At 14 TeV on the other hand, effects are at the 1% level 
and will be hard to see. Because most strings at the LHC are very massive, collapses will 
be rare. Only the drag effect for large rapidities is of any importance, as can be seen from 
the figure. At higher energies, higher order effects will become important - e.g. at the 
LHC, less than half of the total heavy quark cross section comes from the leading order 
fusion processes included in fig. [| 

We also show some examples of charm photoproduction in 7p collisions. In this case 
the study is further complicated because of the photon structure, which can be classified 
as direct or resolved. In the direct case, the photon interacts as a point particle in the 
hard scattering subprocess. If the photon has fluctuated into a qq pair it can interact as a 
meson and therefore you need to know the parton distribution of the photon. These two 
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Figure 6: Charm cross section for direct and resolved photons with different p± cuts. 

classes of events give rise to fairly different event structures, see fig. |6]. The drag effect 
shifts the hadron spectrum in the direction of the proton beam remnant. 



5 Summary 

We have studied some non-perturbative effects in heavy quark production within the 
framework of the Lund string fragmentation model. We find some agreement with data 
if specific parameters are used. We also apply the model to pp and 7p physics and arrive 
at predictions of bottom- asymmetries at HERA-B and the LHC. 
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